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Abstract. We demonstrate that a local mapping / in a space of bise- 
quences over {0, 1} which conserves the number of nonzero sites can be 
viewed as a deterministic particle system evolving according to a local 
mapping in a space of increasing bisequences over Z. We present an al- 
gorithm for determination of the local mapping in the space of particle 
coordinates corresponding to the local mapping /. 



1 Introduction 

Cellular automata (CA) are dynamical systems characterized by discreteness 
in space and time. In general, they can be viewed as cells in a regular lattice 
updated synchronously according to a local interaction rule, where the state of 
each cell is restricted to a finite set of allowed values. Among many applications of 
CA, models of road traffic flow, first proposed by Nagel and Schreckenberg in 1992 



[12], attracted substantial attention in recent years. Many theoretical aspect of 
the Nagel-Schreckenberg model are still not fully understood, and therefore several 
simplified models have been proposed, including models based on deterministic 
cellular automata |11, ^. 

One of the interesting features of these models is the fact that they can be 
described using two equivalent paradigms: either as one-dimensional cellular au- 
tomata or as systems of interacting particles on one-dimensional lattice. The sim- 
plest example is rule 184, one of the elementary CA rules investigated by Wolfram 



1 15], and later extensively studied in the context of surface growth models JlO], as 
well as in the context of density classification problem ]Q. It is one of the only 
two (symmetric) non-trivial elementary rules conserving the number of active sites 
]Q 0], and, therefore, can be interpreted as a rule governing dynamics of particles 
(cars) . Particles (cars) move to the right if their right neighbor site is empty, and do 
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not move if the right neighbor site is occupied, all of them moving simultaneously 
at each discrete time step. Using terminology of lattice stochastic processes, rule 
184 can be viewed as a discrete-time version of totally asymmetric simple exclusion 
process. 

A general question which can be asked is when a given CA rule can be treated 
as a rule governing motion of particles, and how to determine all rules possessing 
this property. Since the number of particles has to be conserved, it is clear that all 
such rules must conserve the number of nonzero sites. In [Q, ^, CA rules of this 
type and their phenomenology were investigated. 

In this work, we will formalize concepts introduced in We will prove that 
each conservative CA can be associated with a local mapping of particle coordinates, 
and we will demonstrate how such a mapping can be constructed. 

2 Bisequence spaces and their mappings 

The set of definitions given below closely follows terminology used in , with 
minor modifications. 

Let 5 be a countable set with cardinality 5, which will be called a symbol set. 
If g is finite, we will assume that Q = {0,l,...,.g — 1}, otherwise we will often 
assume Q — "Z. 

A bisequence over Q is a function on Z to Q. Let X{Q) denote the set of 
bisequences over Q, i.e., X{Q) — Q^. If a; G X{Q) and i G Z, then x{i) will be often 
denoted by Xi. 

Let n £N,n > 0. An n-block over Q is an ordered set X1X2 . . . Xn, where Xi G G. 
The set of all n-blocks over Q will be denoted by Bn{G)- 

Let / be a mapping of Bn{G) into Q. The set of all such mappings, for a given 
n G N, n > will be denoted by F{Q,n). Mapping / will be frequently called a 
local mapping, or equivalently a cellular automaton rule. 

Let / G F{G, n) and let m G N, m > 0. Corresponding to / and to, we define a 
mapping : B„i+n-i{G) Bm{Q) as follows. Let b = bi. ..b 

(G) 

and let Oi ~ f{bibi+i . . . 6i+„_i), where i = 1,2, . . . , to. Then we define fm{B) — a, 
where a — 0102, ■ ■ ■ ,am- 

Finally, let us define a global mapping foo ■ X{G) X{Q) corresponding to 
a local mapping / G F{Q,n). Let x G X(Q), and let y G X{Q) be defined as 
Vi = f{xi, Xi+i, . . . , Xi+n-i) for aU i G Z. Then we define foc{x) = y. 

For a G Z we can define a mapping foe, a ■ X{G) X{Q), similar to foe- If 
X G X{g), and if y G X{g) such that yi = f{xi^a,Xi+i-a, ■ ■ ■ ,Xi+n-i-a) then, by 
definition, foo,a = y- Note that foo = /oo,o- 

As an example, consider a set of mappings F({0,1},3), which were studied 
in |l5| and called elementary cellular automata rules. There is 2^^ = 256 such 



mappings, and it is customary, following [15|, to assign them code numbers C{f) 
defined as 



1 



C{f)^ E /(xi,X2,a;3)2(2^-+2^-+2"-3). (2.1) 

Example 2.1 For instance, the local mapping with code number 184, often 
referred to as simply rule 184, is defined by 

/(0, 0, 0) = 0, /(0, 0, 1) = 0, /(0, 1, 0) = 0, /(0, 1, 1) = 1 (2.2) 
/(1,0,0) 1, /(1,0,1) = 1, /(1,1,0) - 0, = 1. 
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It is straightforward to show that the above definition can be written in a more 
compact form as 

f{xi,X2, X3) = X2+ min{a;i, 1 - X2} - min{a;2, 1 - X3}. (2.3) 

Having the local mapping, we can compute m-block mappings fm- For example, 
/3(10100) = 101 because /(I, 0, 1) = 1, /(0, 1, 0) = 0, and /(I, 0, 0) = 1. 

It is a common practice to view CA as dynamical systems and to investigate 
trajectories of points in the bisequence space X(Q), where by a trajectory of a 
point X G Q we mean the sequence {f!^{x)}^Q, superscript k denoting multiple 
composition 

/^ = /oo(/oo(---/oo(a;))). (2.4) 

v ' 

k times 

Since CA are often used as model of real physical systems of a finite size, peri- 
odic boundary conditions are frequently employed. Using the formalism introduced 
in this section, periodic boundary conditions can be easily incorporated by assum- 
ing that the bisequence under consideration is periodic with some period L, i.e., 
there exist some L G N such that Vi G Z : Xi+L = Xi. Obviously, all bisequences 
belonging to the trajectory of a periodic bisequence are periodic. 



3 Conservative rules 

As in any other dynamical system, symmetry and conservation laws play an 
important role in cellular automata. Additive invariants in one-dimensional CA 



have been studied by Hattori and Takesue ||8|, |14| . They obtained conditions which 
guarantee the existence of additive conserved quantities, and produced a table 
of additive invariants for Wolfram's elementary CA rules. In this work, we will 
consider simplest possible additive invariants, namely the number of active sites 
("active" meaning non-zero). 

Let us start from a simple example, rule 184 introduced in the previous section. 
Let us consider an i-periodic bisequence x G A'({0,1}) and its image under the 
local mapping with code number 184, as defined by eq. ( ^.3| ), y — foo{x). The 
number of active sites in a single period of x (to be referred to as simply the 
number of act ive sites in x) is equal to J2i=i ^i- The number of active sites in y is, 
using eq. (2.3), 



i=L i—L i—L 



'^Vi ='^f{x^+l,X^+2,X^+3) = '^Xi+2 + 

2—1 2 — 1 i — 1 

i—L i—L 

^min{xi+i, 1 - Xi+2} - ^min{a;i+2, 1 - Xi+3}. 



Since for a periodic lattice the last two sums cancel, we obtain 

i—L i—L 

i=l 1=1 

meaning that the number of active sites is conserved. 
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In general, a local mapping / G F{Q, n) will be called L- conservative if for any 
L-periodic bisequence x € X{Q) the following condition is satisfied: 

i—L i—L 

^ /(Xi, Xi+i, . . . , Xi+n-l) = ^Xi. (3.2) 
i=l i=l 

Using the fact that x is periodic, this condition can be rewritten as 

f{xi,X2, ... ,X„) + f{x2,Xz, . . . ,Xn+i) H h f{xL,Xi, . . . ,Xn-l) 

= xi + X2 + ■ ■ ■ + XL- (3.3) 

If for every L > a mapping / e F{Q, n) is L-conservative, than it is said to be 
conservative. The following theorem Q is helpful in determining if a given rule is 
conservative. 

Theorem 3.1 (Q) If a mapping f E F{Q ,n) is L-conservative for L = 2n~2, 
then it is conservative. 



Pr oof To prove the above result, we shall show that, if L > 2n— 2, any equation 
of (3.3) type, i.e., equation expressing conservation condition for L-periodic con- 
figuration, is a linear combination of three equations which express, respectively, 
conservation conditions for {L — 1)-, (2?! — 3)-, and (2?! — 2)-periodic sequences. 
More precisely, for all L-pcriodic bisequcnces with period L, {xi, X2, ■ . ■ ,xl}, Equa- 
tion (|3.3|) can be written 



(/(a;i,X2, ... ,Xn) + f{x2,X3, . . . ,X„+i) H h f{xL-l,Xi, . . . ,Xn-l)^ 

-(^f{xi,X2, ■ . . , Xn-2, XL-n+1, XL-n+2) 

+ f{x2,X3, . . . , XL-n+s) ^ 1- f iXn-2, • ■ ■ , XL-i) 

+ f{xL^n+l,XL~n+2, ■ ■ ■ ,Xl~i,Xi) + . . . + f{xL_i,Xi, . . . , Xn-2, X L-n+l)^ 

+ (^f{xi,X2, ■ . ■ , Xn-2i XL^n+1, XL-n+2) + f{x2,X3, . . . ^XL-n+s) H 

+ f{Xn--2,XL-n+l, ■ ■ ■ ,Xl^i) + f [x L-n+1, X L^n+2 , ■ ■ ■ ,Xl^1,Xl) + ■ ■ ■ 
+ f{xL-l,XL,Xi,... ,Xn-2) + fixL,Xi, . . . ,Xn^2,XL-n+l)j 

= {Xi+X2'^ h XL-l) - (a^l H h Xn-2 + XL-n+1 H h XL-i) 

+ {X1 H h Xn-2 + XL-n+1 H h X l ) ■ (3.4) 

The above rearrangement of term^ is valid if xl-i = xl, and for almost every 
periodic bisequence, we can choose a "coordinate system" (labeling of sites) in which 
it is true. The only exception is, for even L, the bisequence with period 1010 ... 10. 
Verifying (3.4) is in this case possible, but needs to be done with a different method 
(see below). By induction, relation ( ^.4[ ) shows that the conservation condition for 
L-periodic bisequence is a linear combination of conservation conditions for (2n— 3)- 
and (2n — 2)-periodic bisequences. Using similar rearrangement of terms as in ( |3.4D , 
it can be shown that conservation conditions for (2n — 2)-periodic bisequences imply 
conservation conditions for {2n — 2)-periodic bisequences, concluding the proof. 



'^As pointed out by a referee, this rearrangement of terms can be also understood as a 
representation of an L-cycle in a de Bruijn graph as the sum of an L — 1-cycle, a. 2n — 2-cycle, 
and a, 2n — 3-cycle. 
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Consider now the exception mentioned in the preceding paragraph, the bise- 
quence 1010 . . . 10. Conservation condition for arbitrary L reads 

/(1010...10) + /(0101...01) + --- + /(0101...01) = ^, (3.5) 



if n is even, and 



/(1010...01) + /(0101...10) + --- + /(0101...10) = -, 
„ ' 2 



(3.6) 



if n is odd. That is, 



if n is even, and 



/(1010...10) + /(0101...01) = 1, 



/(1010...01) + /(0101...10) = 1, 



(3.7) 
(3.8) 



if n is odd. In the above two equations L does not occur at aU, so obviously, if eq. 
( |3.5[ ) or ( |3.6[ ) is true for i = 2ri — 2, it is true for any L. □ 

Using this theorem, if one wants to find aU mappings / 6 F{Q,n) for some n 
which are conservative, it is enough to check which of them are 2n — 2-conservative. 
In fact, for many mappings it is possible to show that they are not conservative (and 
thus eliminate them from a list of mappings "suspected" for being conservative) by 
employing some of their elementary properties. 



Theorem 3.2 Let n G N and n > 0. If a mapping f e F{0, n) is conservative, 



then 



(a) /(0,0,... ,0) = 0, 

(b) /(1,1,... ,1) = 1, 

(c) E f{xi,^2, 
xi,X2,... ,a;„£{0,l} 



i.e., f is 1-balanced. 



Proof Part (a) becomes obvious if we consider periodic bisequence consisting 
of all zeros, i.e., containing no active sites. Image of this configuration must also 
contain no active sites, and it is only possible if /(O, 0, . . . , 0) = 0. Proof of (b) is 
identical. 

Part (c) says that among all possible configurations of arguments xi,X2, ■ . ■ ,Xn 
of / (there are 2" such configurations), exactly half results in f {xi, X2, ■ . ■ , Xn) — 1, 
as can be seen, for example, in definition of rule 184 (eq. ^.2|) . 

To prove part (c), consider periodic bisequence t G X({0,1}) of period m — 
2". Now, let us construct a set of m blocks A = {^'^-'-'j^'Li such that b^^^ = 
tj, ij+i, . . . , ij+n-i (superscript (j) denotes here just a consecutive number of the 
block b^^^ in the set A). 

Assume that we can find t such that all blocks in A are different. This means 
that each possible block of length n occurs in A once and only once, and therefore 



E f{xi,X2, ... ,Xn) 

xi,X2,... ,a:„G{0,l} 



i=l 



/(i„i,+i,... ,t,+„_i) ^TVi, (3.9) 



where A^i is the number of I's in a single period of t. 
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On the other hand, t G X({0,1}), which is obtained from t by replacing aU 
zeros by ones and vice versa, must also have the same property as t, i.e., in a set 
A = {b^^^jLi such that b'^^^ = ii,ti+i, . . . ,fi+„_i, all elements are different, hence 

m 

X! f{xi,X2,.-. ,Xn) = ^f{ti,ti+l,--- ,tt+n-l) ^ No, (3.10) 

xi,X2,... ,x„e{Q,l} i=l 



where Nq is the number of O's in a single period of t. Comparing ( p.9| ) and (3.10) 
we obtain Nq = Ni = 2"/2 — 2"^^, exactly as required. 

The only problem left is to show that, indeed, for any n > 0, we can construct 
2"-periodic configuration t such that all blocks occurring in a single period of t are 
different (and therefore, constitute a set of all possible blocks of length n). For 
example, for n = 2, consider t with a period 1100. One can easily check that all 
possible blocks of length 2 occur in a single period of t: 11, 10, 00 and 01. For 
71 = 3, t with a period 11101000 has the same property. Again, one can see that 
blocks occurring in the period of t, 111, 110, 101, 010, 100, 000, 001, and Oil, are 
all possible blocks of length 3. 

For a general n, the required t is equivalent to a hamiltonian cycle in the de 
Bruijn graph Q of dimension n (or an eulerian cycle in the de Bruijn graph of 
dimension n — 1). It has been demonstrated that such a cycle always exists (more 
precisely, for a given n, there exist exactly 2^" of such cycles - see, for example. 



review article [13|). □ 



4 Conservative cellular automata viev^red as deterministic particle 

systems 

Since conservative rules conserve the number of active sites, one can identify 
active sites with "particles" which change position after each application of the rule 
/cx3, but their number does not change. Therefore, instead of describing this system 
of particles in terms of lattice sites being occupied or empty, we can describe it by 
specifying list of coordinates of all particles. 

Definition 4.1 An increasing bisequence over Z is an increasing function on Z 
to Z. Set of all increasing sequences over Z will be denoted by X(Z). If s S X^E) 
and i G Z, then s(i) will be often denoted by Si. Now, let us define a mapping 
(f) : X(%) ^ -''^({0, 1}) as follows: let s G X(Z) be an increasing bisequence and let 
X G X({0,1}) be defined by 

1, if3fcGZ:sfe = i, 
0, otherwise, 

for every i G Z. Then, by definition, 0(s) — x. 

Using the notion of particles, s in the above definition can be understood as 
a list of coordinates of particles on one-dimensional lattice. We require this list to 
be increasing, so no two particles occupy the same position. On the other hand, 
x G ^({0, 1}) is just a list of occupancy numbers: for all sites, Xi is either 1 
(meaning site i is occupied) or (empty site). Transformation takes a list of 
particle coordinates s and returns the corresponding list of occupancy numbers. 
For example, if in some s G X{'£) we have s\ = 7, S2 = 9, S3 = 12, then in the 
corresponding x — (f>{s), sites 7, 9, and 12 are occupied, meaning that xj = 1, 
xg = 1, and X12 — 1. Note that since s is an increasing bisequence, we can 
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immediately conclude that the site between 7 and 9 as well as both sites between 
9 and 12 must be empty, i.e., = and xio = xu = 0. 

Since conservative mappings conserve the number of particles, we would now 
like to find a local mapping which transforms coordinate lists y G X{Z) with the 
same effect as conservative local mappings / e ({0,1}), n) transform occupancy 
lists X e X({0, 1}). 

Theorem 4.2 Let n e N, and let f e F{{0, l},n). There exist m 6 N, a G 
and g G F{'L^m) such that 

(i) for every s G X{'L), goo,a{x) G X{'L), 

(ii) 4>9<^,a = foa4>- 



Note that Theorem [4.2| can be represented as 

X({0,1}) 



X(Z) 



(4.2) 



X(Z) 



X({0,1}). 



Before we present the proof, let us consider an example which will clarify the 



meaning of Theorem 4.2 



Example 4.3 Equation 2.3 defined a local function for cellular automaton rule 
184 as 



f{xi,X2, Xs) = X2+ min{a;i, 1 - X2} ~ min{a;2, 1 - X3}. 



(4.3) 



Corresponding rule g of Theorem 4.2 (constructed by an algorithm to be presented 
later) is defined as: 

g{si,S2) Si + min{s2 - si - 1, 1}. (4.4) 



Consider now a bisequence s G ^(Z), in which si = 1, S2 = 2, and S3 = 4, while 
So < 1, S4 > 7. (this means that among lattice sites i = 1 . . . 8 only i = 1, 2 and 4 
are occupied). Let p — goc{s), which means that pi = g{si^ s^+i), yielding si = 1, 
S2 = 3, and S3 — 5. List of particle coordinates {...1,2,4,...} is transformed by 
goo into {...1,3,5,...} 

Now, consider x G -'^^({0, 1}) such that x = 0(s), and let us find y = foc{x), so 
thatQ yi = f{xi-i, Xi, Xi+i). Both x and y can be represented as 





1 


1 





1 















1 





1 





1 











where the first line corresponds to x, and the second line to y. One can readily see 
that in the first line particles are located at sites 1, 2 and 4, while at second line at 
1,3 and 5, confirming that (j)goc — foc4>- Note that in this example application of 
goo results in a well defined motion of particles: if the site on the right of a given 
particle is empty, it moves to that site, otherwise it stays in the same place. Using 
notation introduced in pJ , this can be written as 



10, 



o 

11. 



(4.5) 



^Strictly speaking, we are using faa,l here, not /oo, mainly due to historical convention 
commonly employed for elementary (3-input) cellular automata. If we wanted to use yi = 
f{xi,Xij^l,XiJ^2), we would have to slightly redefine g, so that 9(51,^2) = si — l+min{s2~si — 1, 1}. 
However, the form used in the example is more intuitive (it represents asymmetric exclusion pro- 
cess with discrete time). 
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(arrow shows where the particle wiU move, circular arrow indicates that the particle 
stays in the same place). 

5 Labeling scheme for lattice sites 



In order to construct function g of Theorem 4.2, we will first introduce trans- 
formation of coordinate space X{Z) to intermediate ("mixed") space Vf, which 
combines both coordinates of particles and information about occupancy of lattice 
sites. This can be done by labeling occupied lattice sites with consecutive integers, 
and empty lattice sites with extra symbol •, so that configurations of particles be- 
come bisequences over Z U {•} (the extra symbol "•" had to be introduced since 
"0" can be a particle label, and therefore cannot be used to denote empty sites). 

Definition 5.1 Let x G X(Z). Let y G X(ZU{»}) be a bisequence constructed 
as follows: yi — j if there exists j G Z such that Xj = i, otherwise yi = •. By 
definition, ^l^{x) = y. Bisequences y as defined above will be called increasing 
labeling bisequences, and the set of all such bisequences will be denoted by ^' = 

V(X(Z)). 

Note that ■;/' transforms a bisequence of particle coordinates, such as a; = 
...1,2,4,... into a bisequence ipi^)'- 



which resembles (j){x) , except that occupied sites are now labeled with unique inte- 
gers, increasing from left to right (the label of an occupied site is always larger by 1 
than the label of the closest occupied site on the left). If we define, for G ZU {•}, 

then, obviously, [<f){x))i — \\{'4'{x))i\\ for every i G Z. 

We will now define a mapping which plays a similar role in 5* as / plays in 
X({0, 1}). Corresponding to /, mapping / G F(Z U {•},?!) will be defined as 
follows. 

Define 

G(xi,X2,... ,a:„) = /(0,0,... ,0,||xi||) + /(0,0,... ,||a;i||,||x2||) + 

. . . + /(O, 11x211, . . . , ||x„||) + /(llxill, 11x211, . . . , ||x„||), (5.2) 



and 



H{xi,X2, • . . , Xn, k) — min{i G N : ^ ||a;j|| — k}. (5.3) 
Now, function / is defined as 

/(X1,X2,...,X„) = ( if/(ll-llUI^2||,...,||x„||) = l, 

^ ' ' ' ' \ otherwise, ^ ' 

where I = H{xi,X2, ■ ■ ■ ,Xn,G{xi,X2, ■ ■ ■ ,x„)). 

Example 5.2 For rule 184, as defined in (pT^), applying the above definition 
we obtain for every xi, X2, x^ G Z, 

/(•,•,•) = •^{•.•,Xl) = •,!{•, ^ •,/{•, Xi,X2) = Xi, 

/(xi, •, •) = Xi, /(Xi, •, X2) = Xi,f{xi,X2, •) = •, fixi,X2,X3) = X2. 
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An immediate consequences of the definition of / is X2, . . . = 

/(||xi, II, ||x2, II, . . . , ||a;„, II), and therefore: 

Lemma 5.3 Let n ^ N, n > 0, and let f ^ F{{0,l},n). Then for every x ^ , 

\\foo{x)\\^ f^{\\x\\). 

A crucial property of / is that an image of an increasing labehng bisequence 
under f^o is also an increasing labeling bisequence: 

Theorem 5.4 Let n e N, n > Q, and let f e F({0, l},n). Lf x <E , then 

foo{x) e 

Proof Let us consider an increasing labeling bisequence x and its image y = 
foo{x). Consider two sites in y, yi and yj, i < j, such that all sites between them 
are empty, i.e. yk = * for all i < k < j. We want to show that yj = yi + 1, 
since this would demonstrate that y is indeed an increasing labeling bisequence. 
Obviously, yi = f{xi,Xi+i,... ,a;j+„_i) and yj = f{xj,Xj+i,... ,Xj+n-i)- There 
must be some p, q such that yi = Xp and yj = Xq, which, according to the definition 
of /, must satisfy 

p=i-l + H{xi,Xi+i, . . . ,Xi+n-l,Gix^,Xi+i, . . . ,Xi+n-l)), (5.5) 
q=j -1 + H{Xj,Xj + i, . . . ,Xj+n-l,G{Xj,Xj + i, . . . ,Xj+n-l))- 

As a first step of the proof, we will find relationship between G{xi, Xi^i, . . . , Xi-^-n-i)) 
and G{xj,Xj+i, . . . ,Xj+n-i)). 

Let us denote the number of nonzero sites in the block Xi,Xi+i, . . . ,Xj^i by 
N, i.e., N — \\xi\\ + \\xi+i II + , . . . , +||a::j_i ||. Consider now a bisequence t G 5* such 
that ti = xi for i<l<j + n — 1, and ti = • otherwise, and another one, w € 5*, 
such that ui = xi for j <l<j + n— 1, and ti = • otherwise. Due to the fact that 
/ is conservative, we have 

oo oo 

J2 ii/>/,t/+i,... ,t/+„-i)ii-^+ E 

,Ui+i, . . . ,U;+„_i)||. (5.6) 

/— — oo /— — oo 



Using definitions of t and u, and the fact that /(0,0, ... ,0) = 0, equation (5.6) 
becomes 

||/(», . . . ,»,x,)\\ + ||/(», . . . ,»,a;^,a;^+i)|| + . . . + \\f{xi,Xi+i, . . . ,a:i+„_i)|| + 

\\f{xi+i,Xi+2, ■ ■ ■ ,X.,+n)\\ + ■■■ + \\f(Xj,Xj+i, . . . ,Xj+n-l)\\ + 
\\f{Xj+i,Xj+2, ■ ■ ■ ,»)|| + . . . + ||/(.Tj+„_i,», . . . ,»)|| = N + 

||/(»,... ,»,Xj)\\ + ||/(»,... ,»,Xj,Xj+i)\\ + ... + \\f{xj,Xj+i,... ,a;j+„_i)|| + 

\\f{Xj + i,Xj+2,... ,»)|| + ...+ \\f{Xj+n-l,»,--- ,»)||- (5.7) 

However, we know that there are no particles between yi and yj, therefore 

\\f{xi,xi+i, ... , a;;+„_i)|| = (5.8) 

for all I such that i < I < j . Taking this into account, and using definition of G 
(cq. |3.2[ ) we obtain, after some cancellations, 

G{Xi,Xi+i, . . . ,Xi+n-l) + \\fiXj,Xj+i, . . . ,X-,+„_i)|| = 

N + G{xj,Xj+i,... ,Xj+ri-i)- (5.9) 
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Since ||/(a;j, x^+i, . . . , a::j_|_„_i)|| = 1, this finally becomes 

G{Xi,Xi+i, . . . ,Xi+n^i) = N -1 + G{Xj,Xj+i, . . . ,X.j+n-l)- (5.10) 

Having this relationship, note that, according to definition of H (eq. |5.3| ), 
Xp is fcp-th particle in the block Xi, Xi+i, . . . , Xi+n-i starting from the left, where 
kp — G{xi, Xi+i, . . . ,Xi+„_i). Obviously, it is also is fcp-th particle in the block 
Xi,Xi+i, . . . , Xj+n-i starting from the left. 

Similarly, Xq is fcg-th particle in the block Xj,Xj^i, . . . ,a;j+„_i, where fcg = 
G{xj, Xj+i, . . . , Xj+n-i)- Knowing that sites Xi, . . . , Xj-i contain exactly N parti- 
cles, we conclude that Xq must also be (A^-l-fcg)-th particle in the block Xi, Xi+i, . . . , 
Xj+n-i- However, fcg = fcp — + 1, therefore: 

Particle Xp is kp-th particle from the left in the block Xi, Xi^i, . . . , Xj+n-i, while Xq 
is {kp -\- l)-th particle in the same block. 

This means that there is no other particle between Xp and Xq, hence Xq = .Tp + 1 
(we assumed that x is an increasing labeling bisequence), and, finally, yj + 1, 
which is exactly what what we wanted to show. □ 



6 Construction of local mapping in X(Z) 



We are now ready to see how g of Theorem 4.2 can be constructed. 

First, let us extend definition of the mapping ip (Definition 5.1) to blocks of 
particle coordinates. Let b e BmC^) be called an increasing block if 6i < 62 < . . . < 
bm- Set of all such blocks will be denoted by BmC^)- For m > 1, b G BmC^) let us 
define c — ^{b) such that c G Bb^-bi+i{'^^ {•}) and that, for every I £ N satisfying 
1 < ? < 6,„ — 5i + 1, if there exist fc G N such that bk = I, then c/ = fc — &i -I- 1, 
otherwise q = •. Mapping tp, as before, transforms finite and increasing list of 
particle coordinates into finite block in particles are located at sites bi — bi + 1 and 
labeled with labels I. For example, if 6 = {9, 10, 14}, then ip{b) — {1, 2, •, •, 3}. 

Similarly as we defined increasing labeling bisequences, we define blocks with 
increasing labels as elements of the set ip{Bm{'^))- This set will be denoted by 
i3™(ZU{.}) 



In Theorem 5.4, we proved that foo maps increasing labeling bisequences to 



increasing labeling bisequences. The same is true for blocks with increasing labels 



(and can be proved using almost identical reasoning as in the proof of Theorem 5.4, 
therefore we omit the proof): 

Corollary 6.1 For any m,n E N , m > n, and for f E F{{0,l},n), if b E 
S™(ZU{.}), thenf^{b)EB 

'm— n+1 

(zu {•}). 

Let us now assume that we have a list of particle coordinates s E X{Z), and 
we want to find new list of particles' positions t such that (p{t) — /oo(0(s)), where 
/ E -^({0, l},^^) for some n E N, n > 1. Using /, this can be done as follows: 

Algorithm 6.2 For each particle coordinate Si perform the following steps: 

1. construct a block b = {si_„+i, Si_„+2, . . . , Si+„-i} 

2. find block c = /s;^„_j_s;_„^j+i(^/'(&)), this will be a block with increasing 
labels, with particles in their new positions 

3. find d = 'ip~^{c), this will be a block containing new positions of particles 

4. find new coordinate of particle i, which is simply equal to di 
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Note that the size of the block b is chosen large enough to ensure that particle 
with label i remains in the block after step 2. For many rules, smaller neighborhood 
will suffice. 

The above algorithm provides a way to compute ti given Si and coordinates of 
particles in the neighborhood of Si. This means that it defines a local mapping g 

In fact, it is possible to write explicit expression 
Corresponding to / G F({0,l},n), mapping 



requested by theorem Theorem 4.2 
for g based on the Algorithm 6.2. 
g : /B2n-i(^) ^ Z is defined as follows: 



and parameter a of Theorem 4.2 equals to n 



As a final note, let us remark that the relationship of /, /, and g can be 
represented as 



X(Z) 



IMI 



IMI 



X({0,1}) 



X({0,1}). 



We should also mention that g constructed here is not the only local mapping 
satisfying Theorem 5.4. Infinite number of local mappings with the same properties 
as g can be constructed as a superposition of g and the shift map on X{Z) (by the 
shift map we mean a : X(Z) —>■ X(Z) such that {a{x))i — Xi-i for every x £ X(7j) 
and every i £ Z). All mappings a'^g as well as <T~^g for any k gN could be used in 
place of 5. As a convention, we will always use a local mapping such that if Xi = i 
for every i S Z, and y ~ goo{x)^ then — i, i.e., if all lattice sites are occupied in 
X, then in y they are in the same position. 



7 Example: four-input conservative rules 

Among 2^* local mappings / G -^({0, 1},4) (four-input rules), only 22 rules are 
conservative For all of them, we computed /, as shown in Table 1. 

The first column in this table represents code number C(/), defined similarly 
as (O): 



Cif) 



^ f{x,,X2,xs,x,)2^^'-^^+^'^^+^'--^+'°^^^\ (7.1) 



The second column is a binary representation of C(/), meaning that it is a sequence 
of 16 binary digits 015014 . . . ao such that for every xi,X2, X3, X4 G {0, 1} 



a23xi+22x2+2i2:3+20i4 = f{xi,X2,X3,X4). 



(7.2) 



Entries in the third column are constructed in a similar way. They are sequences 
&15614 ... 60 such that if in the definition of / we have f{xi, X2, X3, 2:4) = Xj, then 
we define 



and if 



^23||xi||+22||a:2||+2i||x3||+20||x4|| - ^j, 
f{xi,X2,X3,X4) = •, 



(7.3) 
(7.4) 
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Table 1 Four-input rules conserving density of nonzero sites. 





Binary form 


/ 


43690 


1010101010101010 


X4 • X4 • X4 • X/i • X4 • X/i • X4 • Xi» 


43944 


1010101110101000 


x^ * X2 * x^ * x\xixz • a;2 • • •• 


47288 


1011100010111000 


Xs • X2X2X'i • • • X3 • 3:2X2X3 • •• 


48268 


1011110010001100 


X3 • X2X2X3X3 • •xa • • • X3X3 • • 


48770 


1011111010000010 




49024 


1011111110000000 




51448 


1100100011111000 


.T3X3 • •X3 • • • X3X3X2X2X3 • •• 


52428 


1100110011001100 


0:3X3 • •X3X3 • •X3X3 • •X3X3 • • 


52930 


1100111011000010 


0:3X3 • •X3X3X4 • X3X3 • • • •X4« 


53184 


1100111111000000 




56528 


1101110011010000 


X2X2 • X2X1X1 • •X2X2 • X2 • • • • 


57580 


1110000011101100 




58082 


1110001011100010 


0:3X30:4 • • • 0:4 • 0:3X3X4 • • • 0:4* 


58336 


1110001111100000 




59946 


1110101000101010 


X3X3X4 • X4 • X4 • • • X4 • X4 • X4« 


60200 


1110101100101000 


X2X2X2 • X3 • XiXi • •X2 • X3 • •• 


61680 


1111000011110000 


X2X2X2X2 • • • •X2X2X2X2 • • • • 


62660 


1111010011000100 


X2X2X2X2 • X3 • •X2X2 • • • X3 • • 


63544 


1111100000111000 




64524 


1111110000001100 




65026 


1111111000000010 




65280 


1111111100000000 





then 



&23||a;i||+22||a;2||+2i|k3||+20||a;4|| 



(7.5) 



This win become clear when we consider, for example, rule 53184 from Table 1. 
For this rule, in the third column we have X2X2 • •X1X1X1X1X2X2 This 
means that for xi, X2, X3, X4 € Z 



/(X1,X2,X3,X4) = • 
/(X1,X2,X3,X4) = • 

/(X1,X2,X3,X4) = X2 
/(X1,X2,X3,X4) = Xi 



if 
if 



|Xi,X2,X3,X4|| 
|X1,X2,X3,X4|| 



if ||X1,X2,X3,X4| 
if ||X1,X2,X3,X4|| 



0000 
0001 

0111 
1000 



/(xi,X2,X3,X4) = X2 if ||xi,. 0:2, X3,. 0:411 = nil, 

which is a set of 16 equations fully defining / for rule 53184. Column 3 of Table 
1 represents sets of similar 16 equations in a condensed form for all conservative 
4-input rules. 

In order to reduce the number of rules for further considerations, we can exploit 
the fact that cellular automata rules obtained by spatial reflection or conjugation 

(interchanging zeros and ones in rule tabic), or both reflection and conjugation 
of a given rule have similar properties as the original rule. Let, for a given / G 



Published in "Hydrodynamic limits and related topics", S. Feng, A. T. Lawniczak 
and S. R. S. Varadhan (Eds.), AMS, Providence,RI (2000), pp. 57-69. 



A class of cellular automata equivalent to deterministic particle systems 



69 



F({0,l},n), 

fRixi,X2,X3,X4) = f{x4,,X3,X2,Xi), (7.6) 

fc{xi,X2,X3,X4) = 1 - /(I - Xi, 1 - X2, 1-X3,1- X4), (7.7) 

fRcixi,X2,X3,X4) = fcix4,X3,X2,Xi). (7.8) 

All rules in the set {/, fji, fc, fnc} have similar dynamics, therefore it is enough to 
consider only one of them. For conservative 4-input rules, we have seven such sets, 
or equivalence classes: {49024, 59946, 65026, 43944}, {53184, 58082, 64524, 47288}, 
{56528, 57580, 62660, 51448}, {58336, 52930, 63544, 48268}, {60200, 48770, 60200, 
48770}, {61680, 52428, 61680, 52428}, and {65280, 43690, 65280, 43690}. We will 
now present mappings g for the first rule in each set, defined in terms of the step 
function O (for to G Z, 8(m) = 1 if m > 0, otherwise 9 (to) = 0) and the delta 
function (for m, n G Z, Sm.n = 1 if to = 1, otherwise Sm,n — 0): 

Rule 49024: g{xo, xi, X2, x^,) — xi — \ + 8(2:3 + X2 — 2xi — 3), 

Rule 53184: g{xQ, xi, X2, x^,) = xi — \ + 8(2:2 — xi — 1), 

Rule 56528: 5(2:0, a:i, 2:2,2:3) = 2:1 - 1 + 52,x2-xi, 

Rule 58336: 5(2:0, 2:1, 2:2, 2:3) = 2:1 — 1 + 8(2:2 — 2:1 — 2), 

Rule 60200: 5(2:0, 2:1, 2:2, 2:3) = xi - 1 + 8(2:2 - 2:1 - 2) - 5i,x2-xiQ{xi - x^ - 1), 
Rule 61680: 5(2:0, 2:1, 2:2, 2:3) = 2:1 — 1, 
Rule 65280: 5(2:0, 2:1, 2:2, 2:3) = 2:1. 

Note that for all the above rules, 5 is a function of only 4 arguments (with the 
exception of rule 60200, it is really a function of 3 arguments xi, 2:2, 2:3, but for the 
sake of uniformity we always used xq, xi, X2, X3). For all these rules, 4>goa,i = (f'foo- 
Note that we could redefine 5's and introduce 5(^0, xi, 2:2, 2:3) = 5(xo, xi, X2, X3) + c, 
where c E Z, and then we would have 4>goo.i = 4'foo.c- 

8 Concluding remarks 

We have demonstrated that a local mapping / in a space of bisequences over 
{0, 1} which conserve the number of nonzero sites can be viewed as a deterministic 
particle system evolving according to a local mapping 5 in a space of increasing 
bisequences over Z. We also presented an algorithm for determination of the local 
mapping in the space of particle coordinates corresponding to mapping /. 

Viewing conservative CA as systems of interacting particles is often very useful 
in solving problems related to CA dynamics, in particular problems which can 
be broadly characterized as forward problems [^: given a CA rule, determine 
(predict) its properties. For instance, it is often easier to characterize the tra- 
jectory {g]x>{^)}i^o of ^ point s G -^(Z), then the trajectory {/^(x)}^g of a point 
X G Ar({0, 1}. An example of such approach can be found in where we consid- 
ered properties of trajectories of configurations in simplified deterministic models 
of road traffic flow. 
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